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ABSTRACT 


~Y 

The  Perron-Frobenius  Theorem  states  that  a  matrix  with 
nonnegative  entries  has  at  least  one  nonnegative  eigenvalue 
of  maximal  absolute  value  and  a  corresponding  eigenvector 
with  nonnegative  components.  this  paper  we  discuss 
generalizations  of  this  celebrated  theorem  that  locate  an 
eigenvalue  of  maximal  absolute  value  and  the  components  of 
a  corresponding  eigenvector  within  a  certain  angle  of  the 
complex  plane  depending  on  the  angle  which  contains  the 
entries  of  the  matrix.  A  complete  description  of  the  2*2 
case  as  well  as  partial  results  for  the  general  case  are  given. 
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1.  Introduction  and  Statement  of  Results. 


The  Perron-Frobenius  Theorem  states  that  a  matrix  with  nonnegative 
elements  has  at  least  one  nonnegative  eigenvalue  of  maximal  absolute  value. 
It  seems  natural  to  ask  for  generalizations  that  locate  an  eigenvalue  of 
maximal  absolute  value*  within  a  certain  angle  of  the  complex  plane  depending 
on  the  angle  which  contains  the  elements  of  the  matrix. 

More  precisely,  let  us  consider  the  family  A^a)  of  all  n  x  n 
complex  matrices  A  =  (a^)  whose  entries  are  contained  in  a  sector 


S(or)  =  {z:  | arg  z  |  S  or;  0  $  or  S  n  fixed}. 


For  each  A  £  A^a)  let  P  =  P(A)  denote  the  minimal  (nonnegative)  angle 
for  which  the  sector  S(0)  contains  an  eigenvalue  of  A  with  maximal 
absolute  value.  Thus,  defining 


0  (or)  =  sup  P(A) , 

A  €  AjjO*) 

we  pose  the  following 

PROBLEM  1.  Give  0fl(a)  as  a  function  of  o  and  n. 
The  Perron-Frobenius  Theorem  states  that 


0Q(O)  -  Op  n  —  1,2,3,...  • 

Another  simple  observation  is  that  for  any  m  <  n  and  A  £  A^flf) , 

we  have  A  ©  0  c  A  (a),  where  0.  is  the  k  x  k  zero  matrix.  Thus, 
n-m  ^  ~n  x 


(1.1) 


1.2 


Obviously, 

(1-2)  P1(a)  =  a; 

so  (1.1)  and  (1.2)  yield 

Pn(«)  *  PjCo)  =  a,  n  =  1,2,3 . 

Since 

C  V01'5  for  >  a’ 

we  also  find  that 


PQ(a)  is  a  nondecreasing  function  of  a. 

In  Section  2  we  evaluate  and  obtain  some  partial  results  for 

arbitrary  n,  as  stated  in  the  following  theorem: 


THEOREM  1. 

(i)  If  n  =  2 


then 


P2(«)  =  < 


a  for  a  g  ^ 

+  J  for 

n  for  o  >  j  . 


(ii)  If  n  2  3  is  odd  then 

Pn(er)  =  n  for  a  >  . 

(iii)  If  n  4  is  even  then 


1.3 


6  (a)  t  a  +  5-i  n  for  5-  <  a  £ 
n  n  -  Za 


2a 


and 


pn(a)  =  n  for  a  >  ^rj  ■ 

Note  the  discontinuity  of  p2(a)  at  a  =  n/**- 

We  also  note  that  part  (i)  of  the  theorem  implies  that  in  general,  A 
in  A^a)  does  not  have  a  square  root  in  A^(a/2). 

Recall  the* ^  the  Perron-Frobenius  Theorem  not  only  proves  that  for 
matrices  with  nonnegative  entries  an  eigenvalue  of  largest  absolute  value 
is  positive,  but  also  that  a  corresponding  eigenvector  can  be  chosen  with 
nosmegative  components . 

With  this  in  mind,  for  each  A  £  A^Ca)  we  denote  y  =  y(A)  to 
be  the  minimal  (nonnegative)  angle  such  that  the  sector  S(y)  contains 
the  components  of  an  eigenvector  corresponding  to  an  eigenvalue  of  maximal 
absolute  value  which  lies  in  S(0(A)).  Hence,  defining 


we  pose: 


y  (a)  =  sup  y(A), 

“  A  €  ^(a) 


PROBLEM  2.  Determine  y_(a)  as  a  function  of  a  and  n. 

-  n  - - -  - 

The  Perron-Frobenius  Theorem  tells  us  that 


Yq(0)  =  0,  n  =  1,2,3,..., 

and  we  obviously  have 


YjCor)  =  0. 


1.4 


Also,  as  for  Pn(o),  it  is  not  hard  to  see  that 


Finally,  since  any  eigenvector  (Xj,...,xn)'  €  Cn  (prime  denoting 
the  transpose)  can  be  rotated  so  that  its  components  are  embedded  in  the 
sector  S(n-7t/n)  it  follows  that 

Yn(a)  3  Sli 

In  analogy  with  the  Perron-Frobenius  Theorem  it  is  perhaps  natural 
to  ask  whether 

Yn(®)  =  Pn(a) »  d  =  2,3,4,..., 

or  whether  there  is  any  nontrivial  bound  for  YQ(ot)  when  a  >  0.  This 
turns  out  to  be  incorrect,  at  least  for  all  n  £  4,  as  statrd  in  the 
following  theorem  which  is  proven  in  Section  2: 

THEOREM  2. 

(i)  If  n  =  2  then 

for  0  <  fl  S  j 
for  a  >  5  . 


1.5 


Y3(°0  £  2a  +  t|  for  0  <  a 

and 


For  n  t  2  note  the  discontinuity  of  Yn(°0  at  a  ~  0. 

We  also  note  that  in  fact,  for  n  =  2,  a  t  n/4;  n  =  3,  at  n/ 12; 
and  n  t  4,  a  >  0;  Ya(°0  attains  its  maximal  possible  value. 

Theorems  1  and  2  and  the  inequality  in  (1.3)  can  be  summarized  in 
the  following  two  tables: 


Table  1.  Cases  n  =  2,3 


W|N» 


n 

n  §  4  even 

n  £  4  even 

n  £  4  even 

n  5  5  odd 

n  £  5  odd 

a 

0<«^ 

n  .  „  <  n 

2n  0  ~  2n-2 

CM 

1 

K  d 

CM 

A 

O 

0  <  a  ^  5_ 

"  2n 

B 

*n<«) 

i 

£  a 

£  a  +  —  n 
n 

71 

£  a 

n 

Y„<«0 

n 

n 

S=i  R 

n 

5ll„ 

n 

Table  2.  The  case  n  £  4 


We  have  not  settled  the  following  questions: 

QUESTION  1.  If  0  £  a  *  n/2n  and  n  £  3,  is  Pn(a)  =  a ? 

QUESTION  2.  What  is  pn(or)  if  n  >  4  is  even  and  n/2n  <  a  %  n/(2n-2)? 
QUESTION  3.  If  0  <  a  5  n/12,  what  is  y3(«)? 


2.1 


2.  Proofs . 

Prc-oi.  of  Theorem  1.  Since  for  all  n  and  a  we  have 

a  5  PQ(cf)  S  n, 

and  since  Pn(<*)  is  a  nondecreasing  function  of  a,  it  remains  to 
obtain  the  upper  bounds 

(2-D  P2(c0  S  a,  a  $  J  , 

(2.2)  p2(a)  Sct  +  5, 

and  the  lower  bounds 

(2.3)  Pn(«)  §a  +  Slln,  n  *  2  even,  ^  <  a  *  ■£—  , 

(2.4)  PQ(a)  £  n,  n  £  3  odd,  a  >  , 

(2.5)  Pn(a)  £  n,  n  £  4  even,  a  >  ^2  * 

We  start  with  the  upper  bounds. 

Since  multiplication  of  a  matrix  by  a  constant  multiplies  the 
eigenvalues  by  the  same  constant,  we  may  assume  that  the  spectral  radius 
of  A,  i.e.,  the  maximal  absolute  value  of  the  eigenvalues,  is  1. 
Moreover,  since  the  eigenvalues  of  the  conjugate  matrix  A  =  (a^ )  are 
the  complex  conjugates  of  the  eigenvalues  of  A,  we  can  assume  that 

A  =  .W) 


is  an  eigenvalue  of  A. 


With  these  assumptions  we  first  prove  that 

(2.6)  0(A)  <  a  +  n/2  for  A  €  A^a),  0  S  a  g  n/2. 

Let  A  =  (ajjj)  €  AjO*)  have  e*8envalues  ^  =  and  X'. 

X  +  X'=  tr  A  =  axl  +  &21  6  S(of) 

and  since 

X'  6  U 

where  U  is  the  unit  disc  centered  at  the  origin,  then 

X*  e  u  n  (-x  +  s(ct)). 

Thus, 

(2.7)  X  +  X'  €  (X  +  U)  n  S(a)  s  T 

so  that  T  is  the  intersection  of  the  unit  disc  centered  at  X  and 
the  sector  S(a).  Therefore,  if 

0(A)  t  a  +  n/2, 

then  T  =  {0},  and  by  (2.7),  X'  =  -X.  Hence,  I  X* !  =1  with 
larg  X' I  <  arg  X  =  0(A)  in  contradiction  to  the  definition  of  0(A) 
Thus  we  have  (2.6),  and  (2.2)  follows. 

Now,  fix  a  with 

0  Si  a  £  n/4. 


and  let  us  prove  (2.1). 


2.3 


and 


Suppose  that  for  some  A  =  (a^)  €  A2(a)  with  eigenvalues 
A' ,  I  A'  I  ^  1,  we  have 


X  = 


(2.8) 


Since  p(A)  <  a  +  n/2, 

then  if  max.  arg  a..  <  a, 

J  B  JJ 


a  <  p(A)  <  a  +  n/2. 

then  T  #  {0}.  Further,  since  P(A)  >  a, 
there  exists  a  positive  number  r  so  that 


max  arg(a . .  a.)  =  a. 

j=l,2  JJ 


We  can  thus  replace  A  by  the  matrix 


A0  =  ife  (A  +  rXI)  e  A2(a> 


whose  eigenvalues  are  A  and  A"  =  (Ar  +  rA)/(l  +  r).  Hence,  |A"|  <  1, 
unless  A'  =  A  =  ^  tr  A  €  S(A)  which  contradicts  the  assumption  that 
P(A)  >  a.  So,  we  may  use  Aq  instead  of  A,  drop  the  subscript,  and 
assume  without  loss  of  generality  that 

(2.9)  arg  an  =  a. 

By  (2.7)  we  have 

(2.10)  an  +  a22  =  tr  A  6  T; 
so 

we  have  ajj  €  T;  so  if 
in  T,  then  (see  Figure  1) 


Also,  by  (2.9),  (2.10),  and  since  a22  €S(a), 
we  denote  by  a^  the  point  symmetric  to  a^ 


2.4 


a'u  =  2el0(cos(P  -  a)  -  an  €  -an  +  T. 
We  may  not  have 

a22  =  air 


For  in  that  case 


A'  =  tr  A  •  A  =  au  +  aj,  -  \  = 


L(2cr-p). 


*11 


hence  U'  I  =  1  and  by  (2.8),  !  arg  A'  I  <  larg  A.  1  which  contradicts  the 

definition  of  0.  Therefore, 


a22  *  all; 


arg  (A  -  a22)  <  arg(A  -  a^), 

and 

(2.11)  arg(A  -  a.1)(A  -  a22)  =  arg(A  -  an)  +  arg(A  -  a22) 

<  arg(A  -  an)  +  arg(A  -  a^)  =  (B  +  a)  +  (n  -  8  +  a)  =  it  +  2a. 

Since  A  €  A  (a)  we  also  have 
~n 

(2.12)  arg(\  -  an)(A  -  a22)  =  arg(A  -  an)  +  arg(A  -  a22)  ^  2p  >  2a, 
and 

(2.13)  “2a  S  arg  a12an  £  2a. 

But  now,  since  0  £  a  £  fl/4,  equations  (2. 11)- (2. 13)  contradict 
the  characteristic  equation  of  A, 


(A  -  a11)(A  -  a22)  =  <j2a21* 


Thus,  assumption  (2.8)  fails;  so  in  view  of  (2.6)  we  must  have 


0(A)  i  a,  0  S  a  <  n/U, 


and  (2.1)  follows. 


Figure  1. 


2. 


The  main  tools  in  obtaining  the  lower  bounds  in  (2. 3) -(2. 5) 
n  x  n  circulants 


-A 


(2.14) 


A  = 


a  .  aA 
n-1  0 


al  a2 


a  »  a 
n-2  n 


an-3  an-2 


i-l  a0  / 


whose  eigenvalues  are 


where 


A.  =  aA  +  a.ui.  +  a„u)?  +  ...  +  a  .U)®"1,  i  =  l,...,r 
j  0  1  j  2  j  n-1  j  *  J  *  ’ 


tt .  =  e2nii/n,  j  =  l,...,n, 
J 


are  the  n-th  roots  of  unity. 


Assume  now  that 


a  >  n/2n, 


and  set  the  a.  in  (2.14)  to  zero  except  for  the  entries 
J 

*B  B  ’  B+l  re  » 


r  >  0,  0  <  e  <  n/2n,  0  £  a  £  n-2, 


where  e  is  chosen  sufficiently  small  so  that  A  €  A^a)  and  i 
determined  later.  The  eigenvalues  are 

\j  =  (am  +  “jVl)u,j»  j  *  1 . .  ' 

and  we  see  that  their  maximal  absolute  value  is  obtained  for  the 


are  the 


will  be 


j  that 


minimizes 


2.8 


larg  am  -  •rgGya^j)  |  =  |arg  an  -  arg  uy-  arg  a^i  = 


=  (3  -  M 

n  a 


Since  0  <  e  <  7t/2n,  this  occurs  at  j  =  1;  thus, 

1^1 1  ^  I  A  j I  ,  j  =  2,... ,n; 

and  consequently 

P(A)  =  |arg  Aj|. 

Moreover,  since  the  argument  of  a  +  u>,a  goes  from  arg  a  = 

n  l  l  n 

arg  w^a^j  =  a  +  n/n  -  e  as  r  goes  from  zero  to  infinity,  then 
arg  Aj  =  argCa^  +  Vm+lK  =  arg(an  +  ui^)  +  2m/a 
runs  through  the  interval 


(2.15) 


Now,  if  n  is  odd  and 


we  pick 


n/2n  <  a  <  n/n, 


so  the  interval  I  in  (2.15)  contains  n.  Hence,  there  exists  a 
A  e  A  (a)  such  that 

P(A)  =  | arg  Xjl  =  n. 


to 


circulant 


</:***  ‘  .  v  ~ 


2.9 


Therefore, 

Pn(<0  £  n,  n  £  3  odd,  ^  <  a  <  £  , 


and  by  the  monotoaicity  of  Pn(«)  in  (1.4)  we  get  (2.4). 
If  a  is  even  and 

n/2n  <  or  £  n/(2n-2), 


we  choose 


so  that  the  interval  I  contains  the  point  a  +  it  -  7l/n  -  2e.  Thus, 

~BD 

we  can  find  a  circulant  A  e  A  (a)  with 

~n 

P(A)  =  |arg  Ajl  =  a  +  n  -  n/a  -  2e; 
and  since  e  was  arbitrarily  small,  we  have 

pn(a)  Scr  +  Slirt,  BJ2  even,  <  or  S  ^  , 

so  (2.3)  is  established. 

Finally,  if  n  is  even  and 

a  >  n/(2n-2), 

then 

a  >  n/2l 


where  £  s  n-1  is  odd.  Thus,  by  (2.4), 


p£(a)  £  n; 


and  since  0n(a)  is  a  nondecreasing  function  of  n,  we  obtain  (2.5).  □ 

Proof  of  Theorem  2.  Since  for  all  n  and  a  we  have 

*  T  " 

and  since  Yn(a)  is  a  nondecreasing  function  of  a,  it  suffices  to  obtain 
the  upper  bound 

(2.16)  Y2(a)  ^  ,  0  <  a  £  J  , 

and  the  lower  bounds 

(2.17)  Y2C«)  +  0  <  cr  £  £  , 

(2.18)  Y3(c0  S  2a  +  3  ,  0  <  a  £  ^  , 

(2.19)  Yn(°)  *  n,  n  l  4,  a  >  0. 

Starting  with  the  upper  bound  in  (2.16),  we  let 

a  £  n/4, 

and  take  any  A  =  (a^)  g  A2(a).  As  in  the  previous  proof,  we  nay 
assume  without  loss  of  generality  that  the  eigenvalues  of  A  are  A 
and  A'  with 


2.11 


A  =  eip(A),  I  A' |  £  1. 


Now  let  (Xj^)'  be  an  eigenvector  corresponding  to  A.  If  one  of 
the  cosponents  vanishes,  the  other  may  be  taken  to  be  1;  so  in  this  case 
y(A)  =0.  If  the  eigenvector  does  not  have  a  zero  coordinate  we  may 
assume  that  it  is  of  the  form  (l,x)'  where  x  f  0,  thus 

(2.20)  Y(A)  =  ||arg  x|. 

Ue  have 


all  +  al2x  =  A* 
a2i  +  a22x  =  Ax» 


and  therefore 
(2.21) 


K  "  all  _  a21 

a12  A  ‘  a22 


Since  a  £  n/ 4,  then  by  part  (i)  of  Theorem  1  we  have  0(A)  £  a 
so  the  interval  connecting  the  origin  and  the  point  A  lies  in  S(a), 
and  the  line  T  through  A  perpendicular  to  the  interval  [0,A] 

(see  Figure  2)  must  intersect  the  positive  real  axis.  Now,  if  a^j 
and  a22  are  located  on  the  right  of  f  then 


laH+  a22l  >  2  |  A|  =  2, 
in  contradiction  to  the  fact  that 


laH  +  a22  I  =  lA+A’I  £  lAl  ♦  |A'I  S  2. 
Thus,  say,  a^  is  on  the  (closed)  left  side  of  f,  and 


2.12 


2.14 


For  the  lower  bounds,  our  tools  are  n  x  n  lower  triangular 
matrices  of  the  form 


where  evidently,  the  eigenvector  belonging  to  X  is  given,  up  to  an 
arbitrary  factor,  by 

(2.22)  (1,  x,  x^ . xn  *)'  with  x  =  . 

With  these  matrices,  we  now  make  suitable  choices  for  the  quantities 
X,  X’  and  a  for  different  values  of  n. 

First,  for  n  =  2  and 

0  <  a  S  n/4, 
set 

X  =  ei0,  V  =  X  -  .  =  ffi, 

where  £  >  0  is  small  and  then  6  >  0  is  chosen  so  small  that  |X' {  <  1 
and  A  €  A^a).  The  vector  in  (2.22)  is  now 

a,,)'.  „.'-»***/**■> t 


thus 


2.15 


i 


| 

i 


Y(A)  =  ||arg  x|  =  a  +  £  -  |  , 


and  since  £  is  arbitrarily  small  we  obtain  (2.17). 


For  n  =  3  and 


0  <  a  £  n/12, 


we  make  the  same  choice  o£  X,  X’,  and  a  as  for  n  =  2,  so 
the  eigenvector  in  (2.22)  is 


(1,  x,  x2)',  x 


e-i(2or+n/2-6) 


Therefore , 

y(A)  =  |larg  x2l  =  2a  +  n/2  -  e, 

and  again,  since  £  is  arbitrarily  small,  ve  obtain  (2.1 8). 

For  n  =  4  and  a  >  0  we  set 

X  =  e  *e,  X*  =  cos  £,  a  =  sin  £,  0  <  £  <  a. 

The  corresponding  eigenvector  in  (2.22)  is  (1,  i,  -1,  -i)'»  so 


and  (2.19)  holds  for  n  =  4. 

Finally,  if  n  >  4  and  a  >  0,  our  choice  is 

X  =  1,  X'  =  1  -  euJj ,  a  =  E,  £  >  0,  uij  =  e2ni^n, 

where  again,  if  £  is  sufficiently  small  then  A  g  (o)  and 


IX'I  < 


2.16 


The  eigenvector  in  (2.22)  is  now 

(1,  u>,  U)J,  Wj'1)’ ; 

so 

»(A)  =  El  71  , 

and  the  proof  of  (2.20)  is  complete.  □ 

Note  that  the  lower  bounds  in  (2.17)-(2.19)  were  obtained  for  matrices 
in  An(«)  whose  eigenvalues  are  aU  contained  in  S(a).  In  fact,  for 
n  g  4  we  found  matrices  A  such  that  y(A)  obtained  the  maximal  possible 
value,  while  all  the  eigenvalues  were  arbitrarily  close  to  1. 


10.  ABSTRACT  (Canffim*  <n  ramn  tf#a  If  Mea»«7  Identity  by  Mock  iwwtafj 

The  Perron-Frobeniua  Theorem  states  that  a  matrix  with  nonnegative  entries  has 
at  least  one  nonnegative  eigenvalue  of  maximal  absolute  value  and  a  correspond¬ 
ing  eigenvector  with  nonnegative  components.  In  this  paper  the  authors  discuss 
generalizations  of  this  celebrated  theorem  that  locate  an  eigenvalue  of  maximal 
absolute  value  and  the  components  of  a  corresponding  eigenvector  within  a  cer¬ 
tain  angle  of  the  complex  plane  depending  on  the  angle  which  contains  the 
entries  of  the  matrix.  A  complete  description  of  the  2x2  case  as  well  as  par¬ 
tial  results  for  the  general  case  are  given.  _  _ 
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